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Abstract. The level splittings induced by a horizontal magnetic field in a parabolic two-dimensional quan- 
tum dot with spin-orbit interaction are obtained. Characteristic features induced by the spin-orbit coupling 
are the appearance of zero-field gaps as well as energy splittings that depend on the electronic state and 
the orientation of the magnetic field in the quantum-dot plane. It is suggested that these quantum-dot 
properties could be used to determine the Rashba and Dresselhaus spin-orbit intensities 

PACS. 73.21.La - 73.21.-b 



1 Introduction 

Nowadays, spin-related physics has become one of the 
most active branches of research in condensed matter. 
The fundamental physics involved and the potential ap- 
plicability in semiconductor device technology constitute 
the main reasons encouraging this research. In particular, 
there is a special interest in the properties of electronic 
spins confined in a quantum dot. The reduced dimension- 
ality of quantum dots makes them good candidates for sys- 
tems sustaining long-lived spin states ^ and, ultimately, 
allowing coherent spin manipulation ^ . 

A most convenient way to distinguish spin states in 
a quantum dot fabricated within a two-dimensional elec- 
tron gas (2DEG) is by means of a magnetic field. When 
the field is oriented parallel to the plane of the 2DEG it 
is expected that up and down spin directions will split 
by the Zeeman energy, with no additional modification of 
the electronic states. This assumption relies on the com- 
plete decoupling of the orbital motion and the parallel 
magnetic field. Since orbital motion is restricted to the 
2DEG's plane, the vertical Lorentz force associated with 
a horizontal field becomes irrelevant. However, a qualita- 
tive difference appears when the electron feels an impor- 
tant spin-orbit (SO) interaction during its orbital motion. 
In this case the SO mechanism is effectively coupling the 
electronic states with the horizontal magnetic field, be- 
yond a pure Zeeman splitting, even in two-dimensional 
dots. 

Recently, several groups have measured the level split- 
tings of lateral GaAs quantum dots in a parallel magnetic 
field OEl, even for a single electron occupation of the dot. 
The measured splittings are usually interpreted in terms 



of an effective ^-factor in a simple Zeeman formula. It is 
the aim of this work to analyze the theoretical prediction 
of level splittings for a quantum dot in a parallel magnetic 
field when SO interaction is present, emphasizing the de- 
viations from the Zeeman scenario. It is also our purpose 
to discuss the possible SO signatures in the recent exper- 
imental data for single-electron dots 



2 Quantum dot model 

We consider a 2D representation of the effective mass 
Hamiltonian for the GaAs conduction-band electrons (m* = 
0.067me) 
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The above hypothesis implies that the direction perpen- 
dicular to the 2DEG is strongly quantized in comparison 
with the in-plane degrees of freedom. The potential re- 
sponsible for the confinement in the plane is taken as an 
isotropic parabola of constant wq- We also include the Zee- 
man energy 

Hz = l9*fiBB a, (2) 

where B = {Bx^By) and cr = {ax,'^y) are the magnetic 
field and the Pauli- matrix vectors, respectively; /is is the 
Bohr magneton and g* = —0.44 is the bulk GaAs (^-factor. 

The SO interaction is taken into account by adding the 
linear Dresselhaus and Bychkov-Rashba (Sj terms for 
conduction band electrons in a [001] 2DEG |7], 
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These contributions originate, as a relativistic effect, in 
the electric fields present in the heterostructure. In their 
intrinsic reference frame, the moving electrons feel these 
electric fields as effective magnetic fields that interact with 
their spin. The Dresselhaus term is due to the bulk in- 
version asymmetry of the GaAs crystal and its intensity 
depends on the expectation value of the vertical momen- 
tum and a material-dependent constant (7 = 27.5 eVA'^ 
for GaAs) as 



A, 



TT 
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where zq is the vertical width of the 2DEG. On the other 
hand, the Rashba contribution stems from the asymmetry 
of the heterostructure profile (the built-in electric field) 
and it is also sensitive to the electric fields induced by ex- 
ternal gates [Hj. Below we shall treat Ad and A^ as varying 
parameters and study the results for different values. From 
Shubnikov-de Haas as well as magnetotransport measure- 
ments in 2DEG's values for these parameters ranging from 
Ri 5 meVA to ~ 50 meVA have been inferred in GaAs(2l 

nni. 

Adding all the above contributions, the full Hamilto- 
nian TL can be considered as composed of two parts. The 
first one (7io)i containing the kinetic energy, the confin- 
ing potential and the SO interaction, remains invariant by 
the time-reversal symmetry. Conversely, the second part 
of the Hamiltonian, corresponding to the Zeeman term 
'Hz^ breaks this symmetry. More specifically, we have 



H ^ Ho + Hz , 
Ho — Hxy + Hd + Hr 

3 Analytical approximations 



(6) 
(7) 



The analytical diagonalization of the full Hamiltonian is 
not available, although several analytical results can be 
obtained when different pieces of the Hamiltonian are ne- 
glected or treated within perturbation theory. The per- 
turbative analysis when Zeeman and SO terms are of the 
same order of magnitude is rather involved and this case 
will be considered only numerically. In this section we shall 
consider the two limiting cases of zero Zeeman energy with 
a relatively weak SO interaction, and of a Zeeman energy 
much larger than SO effects. 

Of course, another limit with analytical solution ap- 
pears when the Zeeman energy is included and the SO 
interaction is switched off. Under these assumptions the 
Hamiltonian is trivial and its level structure is that of a 
two-dimensional oscillator 

Snis = (271 + Kl -f l)nuJo + ^9>bBs , (8) 

whose states are spin-split by the Zeeman energy gap 
As = \g*\nBB, independently of the orbital state con- 
sidered. In Eq. n = 0, 1, . . . and £ = 0, ±1, . . . are the 



principal and Lz quantum numbers, respectively, while 
s = ±1 is the spin label. 



3.1 Limit of weal< SO in zero field 

Assuming Hz = a perturbative calculation for Hr^d << 
huJo yields corrections of second order in the A's and, due 
to the oscillator and spin degeneracies, requires degenerate 
perturbation theory. The modified energy levels read 



Enes = (2n + \e\ + l)hLUo 



+ ^r) + ^i>^i - ^R)es 



(9) 



The level structure given by Eq. © corresponds to 
that of a two-dimensional oscillator with a constant energy 
shift and a fine structure depending on £ and s. Both shift 
and splitting arc proportional to the SO intensities. It is 
worth stressing that even in the present case of negligible 
Zeeman energy the SO interaction yields a spin-splitting 



of the energy levels given by As = 2^ | £{X 



A?j) 



This zero-field splitting is consistent with the Kramers 
degeneracy present in a half-integer-spin system with time 
reversal symmetry. Indeed, each eigenstate characterized 
by the values of (n, £, s) has a conjugate {n, —s) at the 
same energy. 

An alternative method to derive Eq. Q was suggested 
in Ref. |11| and uses a unitary transformation of Ho lead- 
ing to a new Hamiltonian that is diagonal in spin space 
up to second order in the SO intensities. The remaining 
non-diagonal terms of 0{\^) are small compared with the 
diagonalized ones for typical GaAs A's, and so they can 
be dropped without significant error. Using this alterna- 
tive transformation we analyzed quantum dot properties 
like the far-infrared absorption in Ref. |12| and the spin 
precession in Ref. ^3 ■ 



3.2 Limit of weal< SO in large field 

We assume now that there is a large horizontal field, such 
that Hr,d « Hz- The unperturbed energy levels are 
those of Eq. where up and down spin is defined in 
the direction of the magnetic field B = _B(cos6', sin0). 
As before, a second-order calculation within degenerate 
perturbation theory yields 
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where el!^^^ is given by Eq. (jSJ, and we have defined z = 
g* (irB /huJo, the ratio of Zeeman energy to external con- 
finement frequency, as well as the two auxiliary quantities 



G 
F 



Xr + Xj) -2XrXd sin (20) , 
Xr + XI + 2XrXd sin (26*) . 
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Fig. 1. Magnetic-field evolution of the first three oscillator 
shells of a parabolic quantum dot with SO interaction. The 
magnetic field is along the direction 9 — 4:5 deg and the SO 
intensities are A_d ~ 0.2hujo£o and Xr — O.lhujolo- The labels 
indicate the values of {n,£, s) for each B — level (see text). 



In Eq. the degeneracy of each major shell, with a 
given N = 2n + \£\, is broken into multiplets of states, as 
will be further clarified below when discussing the numer- 
ical results. It can be shown that the energies for n = 0, 
negative £ and s = -1-1 smoothly converge to the exact 
values for large enough magnetic fields. For other values 
of these labels there are discontinuities in the range of va- 
lidity of Eq. H10|l since an additional restriction is that the 
field, besides of being large enough, should not be close to 
a crossing point with an integer z. Note that Eq. ifTUIl may 
actually diverge for z = 1. Another relevant feature of Eq. 
pOfl is that it takes into account the angular anisotropy 
through the sin (26) contribution to F and G. 

It should be emphasized that although the energy lev- 
els from the above perturbative calculations, Eqs. © and 
p(J|l . can still be classified using the uperturbed labels 
(n, £, s), the corresponding states are no longer eigenstates 
of orbital and spin angular momenta. These results gen- 
eralize those of Ref. JHl by including the two sources of 
SO coupling as well as the arbitrary angular orientation 
of the horizontal magnetic field. 



4 Numerical solutions 

As stated above, when both SO and Zeeman terms are 
included on the same footing a numerical resolution is 
required. We have used two alternative methods: a) by 
spatial discretization in a a uniform two-dimensional grid, 
and b) using an oscillator basis. In both methods we do 
not impose any symmetry restriction in either real or spin 
spaces and we have checked that the above analytical lim- 
its are recovered with high accuracy. 

A natural unit system for the calculations is given by 
the parabola energy hujQ and the oscillator length £q = 
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Fig. 2. Angular anisotropy of the first and second shells for 
a magnetic field of modulus B — 20B and the SO parameters 
displayed in each panel. The angle 9 is the polar angle of the 
magnetic field vector. 



h/{m*uJo). The associated unit of magnetic field is B ~ 
uJom*c/e while the SO intensities are given in units of 
UljJqIq. Assuming, for instance, a GaAs parabolic dot with 
hu){) = 1 meV and A^i = 50 meVA we would then have 
4 ~ 338 A, S « 0.58 T and, thus, Xd « 0.15 ^wo^o- 

Figure 1 shows a typical dependence of a quantum 
dot's level structure with an applied horizontal magnetic 
field for the three lowest harmonic oscillator shells. The 
lowest subband corresponds to states having [n — — 
0, s = ±1) and it can be seen that these two states split 
linearly with the applied magnetic field, in a similar way 
to a usual Zeeman splitting. 

The second subband is composed of four states char- 
acterized by (n = 0,^ = ±l,s = ±1). These states show 
a zero-field spin-splitting and at large fields they group 
into two branches, each of them composed of two almost 
parallel close lines. In the intermediate-field regime the 
evolution of these two branches is nonlinear with B. Note 
that the labels are given in the left part of the diagram 
according to the ordering of Eq. (|5J), which is appropriate 
to the 5 = limit. At high i?'s the ordering, given by Eq. 
(|Tn| . changes since s = +1 (—1) corresponds to the lower 
(upper) branch. 

The third subband includes both linear and non-linear 
splittings since this subband is composed of two states 
[n = 1,1 — ^ ±1) that split linearly with B, and 
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four states {n = 0,i = ±2, s — ±1) with a zero-field spin- 
splitting and a non-linear evolution. In general, at low 
fields states with £ — always split linearly with B (thus 
in a Zeeman-like way) while states having ^ ^ do not 
split to first order in the field. 

It can be easily checked that the the analytical ap- 
proximations discussed in Sec. Ill reproduce the zero-field 
splittings and the fine-structured branches at large B in 
Fig. 1. The representation in this figure is a generaliza- 
tion to the parallel case of the well-known Fock-Darwin 
diagrams for quantum dots in perpendicular fields. When 
calculating this type of diagrams it is necessary to specify 
the orientation of the parallel magnetic field in the 2DEG 
plane since, in general, the anisotropy of the SO terms 
can reflect in an angular dependence of the level struc- 
ture. However, there is one case where the level structure 
is fully isotropic and it is when only one source of SO 
coupling is present (Dresselhaus or Bychkov-Rashba). A 
general proof of this statement can be found in the follow- 
ing. Note, however, that the perturbative result Eq. H10(l 
already predicts that the level structure is ^-independent 
when one of the two A's vanishes. 

If only one source of SO interaction is present in the 
Hamiltonian, its part preserving time reversal symmetry 
(7^0 ) fulfills another continuous symmetry S with gener- 
ator Q; i.e., S — eT^^'^l^ . For the Bychkov-Rashba and 
Dresselhaus cases it is — + Sz and G = — Sz, 
respectively. This result follows inmediately from the fact 
that both Ti.R and Hd commute with their correspond- 
ing symmetry generators. Taking into account that the 
energy levels of H and S'^TCS are identical and that the 
latter corresponds in fact to a rotation of the magnetic 
field around a vertical axis and angle 9 it follows that the 
level structure does not depend on the parallel field orien- 
tation. The transformed Hamiltonian indeed corresponds 
to a rotated field since 

S+HS = no + ^g*fiB B ■ (e*«^- a e"'^^-) 

= ^0 + ^.g*Ms Br (T , (13) 

where Br = [B^ cos 9 + By sin 9, By cos 9 — B^ sin 9) is the 
rotated field. This argument breaks down when both SO 
interaction sources are present since, then, the continu- 
ous symmetry S is lost and as a consequence, the level 
structure is anisotropic. 

In agreement with the above discussion, the numerical 
calculations show that the angular anisotropy of the en- 
ergy levels is maximal when Xn ^ Xd and negligible when 
one SO intensity is much smaller than the other. This is 
clearly seen in Fig. 2, which displays the evolution with 
the magnetic field orientation of the levels corresponding 
to the lowest and the first-excited subbands of a dot in 
a parallel field B = 20B. It is worth mentioning that not 
only the relative weight Xr/Xd ^ 1 yields an important 
anisotropy of the level structure. The abolute values are 
also relevant since, obviously, when both intensities are ex- 
ceedingly smaller than Tlluq^q SO effects become negligible 
and the level structure is isotropic irrespectively of the rel- 
ative weight. Another interesting feature is that the level 



anisotropy is sensitive to the relative sign of Xr and Xd- 
This is at variance with most physical properties, like the 
6'-averaged values of the energy levels, that are depending 
only on the absolute values \Xr\ and \Xd\- Focussing, for 
instance, on the first shell gap, we find from the numeri- 
cal calculations that if AfiAu > the gap minima are at 
= 7r/4 + mTT (with m an integer) while if XjiXd < the 
gap minima are shifted to 9 = 37r/4 -I- tott. In all cases, 
the anisotropy is characterized by a periodicity in 20, in 
agreement with the analytical formula derived above. 

The detailed variation of the first-shell splitting with 
the SO intensities when B = 2QB is summarized in Fig. 
3. We note that for the chosen range of A's the devia- 
tion of the 0-averaged splitting Ag from the Zeeman value 
(0.2948fiwo at this B) ranges from to « 20%. The an- 
gular anisotropy (dashed lines) is measured by the differ- 
ence between maximum and minimum energy gap, i.e., 
[ma3i{As{9)) — imx\{As{9))] and, in agreement with the 
above discussion, reaches maximum values along the di- 
agonal fine |Afl| = I Ad I, ranging from zero to ~ 50% of 
As- Figure 4 shows similar results for the splitting of the 
second-shell. Notice that in this more complex case we de- 
fine two splittings: one gives the energy difference between 
the higher and lower states of the second shell while the 
other is associated with the two intermediate ones. 

The numerical results in Fig. 3 qualitatively agree with 
the prediction from Eq. (|10() of an averaged splitting and 
anisotropy proportional to (A|j -I- A|)) and XrXd, respec- 
tively. Nevertheless, the diflference between the actual val- 
ues from the second-order and numerical calculations in- 
creases with the A's and it leads to sizeable errors for val- 
ues above 0.2. For the second-shell results (Fig. 4) the 
errors from the perturbative calculation are more conspic- 
uous since they lead to qualitatively diflferent contour lines 
for A's above ~ 0.15. 
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Fig. 3. Contour plots with the dependence on the spin-orbit 
intensities of the first shell-splitting at B = 2QB. Solid lines 
show the 0-averaged splitting As in Tiuio units (the Zeeman en- 
ergy in these units is Ez = 0.2948). Dashed lines correspond to 
the amplitude of the 9 oscillation, [max(Zis(6)) — ram{As{9))]. 
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Fig. 4. Same as Fig. 3 for the second-shell splittings. Lower 
panel corresponds to the energy gap between higher and lower 
states, while upper panel shows the gap between the two in- 
termediate ones. 



The above results suggest that a high-precission mea- 
surement of the quantum dot level structure could de- 
termine the SO intensities. One should know the specific 
values of the underlying 2DEG (like g* and m*) as well 
as the value of wq. The zero- field splitting of the second 
shell would fix |A|) — A^|. This condition along_with the 
results of the first-shell splitting at large B {As in Fig. 
3) would yield the lower and greater SO intensities in ab- 
solute value, i.e., |A<| and |A>|. The angular anisotropy 
of the splittings could fix, as mentioned above, the rel- 
ative sign. To discern which source (Bychkov-Rashba or 
Dresselhaus) is the greater or lower SO intensity seems a 
somewhat complicate task. In principle, however, it could 
be accomplished by looking at the second shell splittings 
since when \Xd\ > \Xji\ the upper (lower) states of the 
second shell have parallel (antiparallel) spin and orbital 
angular momenta at zero field. 

The role of the orbital level spacing is quite important 
for understanding the i3-evolution of the spin splittings 
when SO interaction is present. To see this, we display in 
Fig. 5 the lowest subband's splitting when the SO inten- 
sitites are kept constant and the confining frequencies are 
varied. Small dots, characterized by a wide gap between 




B(T) 

Fig. 5. Evolution with B of the lower shell splitting for fixed 
SO intensities and different values of the confinement huio- The 
magnetic field orientation in the 2DEG's plane is given by a 
polar angle of S = 45 deg. The dashed line shows the Zeeman 
energy \g*\fj.BB. 



orbital subbands, show no effect of the SO interaction in 
this splitting, but, as the orbital and spin energy scales 
become comparable, due to the proximity of the orbital 
subbands, SO interaction induces a level repulsion that 
reflects in a compression of the levels reducing the spin 
splitting. For the particular A's and magnetic field orien- 
tation in Fig. 5, the dots characterized by an orbital level 
spacing lower than « 1 meV there is a sizeable reduc- 
tion of the splitting with respect to the Zeeman value and 
thus their effective g-factor will be lower than that of the 
2DEG. 



The B-evolution corresponding to the second shell is 
more complex since, as already discussed above, this sub- 
band involves four different states. Figure 6 shows the two 
energy splittings corresponding to the gap between inter- 
mediate states and from the lowest to the highest one, re- 
spectively, of the second shell. In all cases, the evolution of 
the splittings is clearly non-linear showing a positive cur- 
vature at low B's and the above-mentioned finite value at 
zero field. It can be seen that as the orbital level spacing 
is reduced the values of the splittings at large B are also 
reduced due to the level repulsion induced by SO inter- 
action, similarly to the case of the first subband. A very 
interesting feature is that for low enough Hluq the SO cou- 
pling induces a large separation of the two i?-dependent 
gaps. A difference that does not appear in the case of small 
dots (high wo's) and, thus, it reveals the SO-induced fine 
structure of the dot level spectrum. 
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A^ = 30meVA /lQ = 60meVA 




5 10 5 10 5 10 15 
B(T) B(T) B(T) 



Fig. 6. Same as Fig. 5 for the energy splittings of the second 
shell. One splitting gives the energy gap between the highest 
and lowest states of this shell, while the other corresponds to 
the difference between the two intermediate states. 

5 Discussion on experimental evidences and 
conclusions 

In a recent work Potok et al have measured the level 
structure and the corresponding spin splittings of the low- 
est and first-excited subbands for a GaAs lateral dot in a 
horizontal magnetic field (Fig. 2c of Ref. 0). The data 
show a small deviation of the first-shell splitting at large 

from the pure Zeeman result, with a fitted value of 
\g\ ~ 0.37 (to be compared with 0.44 for the bulk). A 
single value of the second-shell splitting at large By is re- 
ported, since the resolution does not seem enough to dis- 
criminate finer structures. This second-shell energy gap is 
found to be somewhat lower than that of the first-shell. 
Unfortunately, the angular dependence is not discussed 
and the zero-field splitting of the second shell is not re- 
solved, although the data suggest the existence of unre- 
solved structures. The value of the confining frequency in- 
ferred from the level spacing is Ucuq « 0.8 meV. As shown 
in Figs. 5 and 6, SO-induced effects for A's in the range 
30 to 60 meVA can reproduce the observations, but only 
for a selected magnetic field orientation. The correspond- 
ing angular anisotropics inferred from Figs. 3 and 4 are 
important and would reduce/enhance the difference from 
the pure Zeeman value when varying the polar angle 6. 

Hanson et al have also measured the spin-splitting 
of the first shell of a GaAs parabolic dot, this case evidenc- 
ing a clear nonlinear -dependence (see Fig. le of Ref. 
P]). The deviation from linear behavior is qualitatively 
similar to the Tiluq = 0.5 meV results of Fig. 5 although 
the experimental value of Tiloq at low i^y 's is close to 1 
meV. Therefore, to ascribe the observation to SO effects 
one should assume that the value of the confinement tiujo 
decreases as increases [T^. 

In summary, we have analyzed the energy level struc- 
ture and splittings that are predicted within the effective- 
mass model of a 2D parabolic quantum dot with Bychkov- 
Rashba and Dresselhaus SO interactions. The deviations 
from the simple Zeeman-splitting scenario have been stressed 
Characteristic features of the SO interactions are the ap- 
pearance of zero-field splittings, effective g-factors smaller 



than the bare one and sizeable anisotropics when both SO 
sources are of the same order. It has been suggested that 
these features could be used to determine the SO intensi- 
ties from quantum dot measurements. 

Spin-orbit intensities in the range of « 50 meVA yield 
effects that are close to recent observations in GaAs quan- 
tum dots. However, a conclusive evidence allowing to val- 
idate or discard the relevance of the SO interaction for 
these measurements is lacking. In this respect, the clarifi- 
cation of the experimental zero-field splitting of the second 
shell as well as of the 6'-dependence of the splittings would 
be very important. 

This work was supported by Grant No. BFM2002- 
03241 from DGI (Spain). 
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